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^ALJSTRACT: The present work involves the discrete modelling of corrosion induced cracking and its influence 

cdbn the bond between reinforcement and concrete. A lattice approach is used to describe the mechanical inter- 
section of a corroding reinforcement bar, the surrounding concrete and the interface between steel reinforcement 
^ind concrete. The cross-section of the ribbed reinforcement bar is taken to be circular, assuming that the inter- 
action of the ribs of the deformed reinforcement bar and the surrounding concrete is included in a cap-plasticity 
f^nterface model. The expansion of the corrosion product is represented by an eigenstrain in the lattice elements 

Cforming the interface. The lattice modelling approach is applied to the analysis of corrosion induced cracking 
i and its influence of the bond strength. The model capabilities are assessed by comparing results of analyses 
"yvith those from unconfined pull-out tests reported in the literature. Future work will investigate the influence 

g)f the stiffness of interface elements and the effect of lateral confinement on corrosion induced cracking. 

c^Cey words: lattice, concrete, cracking, reinforcement, corrosion 



-'i INTRODUCTION 



Tie present work involves the modelling of corrosion 
Qnduced cracking of reinforced concrete by means of 
a three-dimensional lattice approach. Corrosion of re- 
inforcement involves the transformation of steel into 
^jjpxpansive rust. If the expansion of rust is restrained, it 

• results in radial pressure in the confining material. For 
reinforced concrete, this radial pressure and accom- 
panying transverse tensile stresses may cause crack- 
ing (lAndrade et al. 19931) . This cracking is not desir- 
^^ble, since it reduces the anchorage capacity of the re- 
inforcement (ILee et al. 20 02). Most of the anchorage 

^Capacity of deformed reinforcement is provided by 

• ^ibs on the surface of the bar, which resist the slip be- 
tween concrete and reinforcement by transferring in- 
clined radial forces into the concrete ( [Tepfers 1979| ). 
The capacity of the concrete to resist these forces 
can be significantly reduced by the corrosion induced 
cracking. Consequently, there is a considerable inter- 
est in developing models which can predict the mech- 
anism of corrosion induced cracking and its influence 
on the bond capacity of reinforced concrete. 

Many of the present models include the effect of 
corrosion induced cracking on bond by reducing the 
bond strength of the interface between concrete and 
steel (ILee et al. 20021) . In these models, the relation- 
ship between the amount of rust and the reduction of 



bond strength is determined empirically. Thus, these 
models are of limited validity for the prediction of the 
influence of corrosion on bond. 

Only very few models describe the expansion of 
the rust, the radial pressure and the transverse stresses 
on the concrete explicitly ( [Lundgren 2005b| ). These 
models have the potential to establish an analytical 
relationship between the expansion of rust, crack- 
ing and spalling. They can be combined with realis- 
tic bond models ( [Lundgren 2005 a| ), so that the influ- 
ence of corrosion induced cracking on the bond ca- 
pacity can be predicted without the need of empirical 
relationships. However, this modelling framework is 
computationally demanding, since it requires three- 
dimensional modelling of the mechanical response of 
the concrete, the bond between reinforcement bar and 
concrete, and the reinforcement bar itself as shown in 
Figured] 

In earlier work, three-dimensional continuum 
based models were used to describe the cracking 
of concrete surrounding the reinforcement. However, 
three-dimensional continuum modelling of cracking 
in concrete is challenging, since it is not straight- 
forward to include the localised deformations into 
the continuum description. Combined with the mod- 
elling of the bond between concrete and reinforce- 
ment, it can be exceedingly difficult. This might ex- 
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Figure 1 : Three-dimensional modelling of reinforced con- 
crete: Concrete, reinforcement and bond between concrete 
and reinforcement are considered as individual phases. 
Left: Concrete cube (light grey) with reinforcement bar 
(dark grey); right: detail of interface between reinforce- 
ment and concerte. 

plain the small number of models which are available 
based on this three-dimensional approach. Discrete 
approaches, such as lattice and particle models, might 
be a favourable alternative to describe this discontin- 
uous problem. Their potential to model corrosion in- 
duced cracking and its influence on bond is assessed 
in the present study. 

Lattice approaches have been used successfully 
in the past to model the failure of concrete, 
as reported by [Sc hlange nand van Mier (1992)] and 



B oland er and Saito (1998)[ The model by Bolander 
has been shown to accurately reproduce analytical 
solutions for elasticity and potential flow problems 
( |Yip et al. 2005} IB olander and Berton 20041) . Further- 
more, it allows for the use of constitutive models, 
which are formulated by means of tractions and dis- 
placement jumps, as commonl y used in interface ap - 
proaches for concrete fracture (ICaballero et al. 2006b . 
These have been shown to result in an element- size 
independent description of crack-openings. This lat- 
tice approach is used in the present study to describe 
the mechanical response of three phases, namely rein- 
forcement, concrete and bond between reinforcement 
and concrete. In this approach, the lattice elements 
do not represent the meso- structure of the materi- 
als (IZubelewicz and Bazant 1987b . Instead, they are 
used to discretise the continuum. Thus, constitutive 
models are required for all three phases, which are 
described in more detail in the following sections. 

2 MODELLING APPROACH 

The present lattice approach for the modelling of cor- 
rosion induced cracking follows the framework de- 
veloped by Bolander and his co-workers. The nodes 
of the lattice are randomly located in the domain 
to be analysed, subject to the constraint of a min- 
imum distance (IZubelewicz and Bazant 1987b . The 
arrangement of the lattice elements is determined 
from the edges of the tetrahedra of the Delaunay 



tessellation based on the randomly placed nodes. 
The cross-sectional properties of these elements are 
obtained from the dual Voronoi tessellation of the 
same set of random nodes. For the interface be- 
tween reinforcement and concrete, the lattice nodes 
are not placed randomly but at special locations, 
such that the middle cross-sections of the lattice el- 
ements form the boundaries between the two phases 
(IBolander and Berton 20041) . The nodes of the lattice 
elements have six degrees of freedom, namely three 
translations and three rotations. These degrees of free- 
dom are related to three displacement and three ro- 
tation discontinuities at the centroid of the middle 
cross-section of the elements. The three rotation dis- 
continuities are related to moments by elastic relation- 
ships. The three displacement discontinuities are used 
in interface constitutive models to determine the cor- 
responding tractions. In the present study, an elastic 
interface model is used for the reinforcement. One 
limitation of the present lattice approach is that it 
only predicts Poisson's ratios of less than 1/4 for 
3D and less than 1/3 for 2D. This is restrictive for 
the 3D modelling of the steel reinforcement, which 
has a Poisson's ratio greater than 1/4. This limitation 
could be overcome by combining the lattice model 
with continu um tetrahedra as discussed for the 2D 
case by |Grassl et al. (2006)] However, this approach 
is beyond the scope of the present study. The inter- 
face model for concrete is based on a combination 
of plasticity and damage, which describes the soften- 
ing and reduction of the unloading stiffness in ten- 
sion as well as the nonline ar hardening response in 
high confined compression (IGrassl 20091) . For the in- 
terface between concrete and reinforcement, a new 
cap-plasticity model is proposed which is described 
in more detail in the next section. 

2.1 Elasto-plastic cap model for the bond between 

concrete and reinforcement 
In the lattice model, the nodal degrees of freedom are 
related to displacement jumps at the middle cross- 
section of the lattice element. This three-dimensional 
displacement jump u c = (u n ,u s ,u t ) T is transformed 
into strains e = (s n , e s , s t ) T by means of the interface 
thickness h as 

h 

The three subscripts n, s and t denote the normal and 
two tangential directions in the local coordinate sys- 
tem of the interface (Figure Ek). The thickness of the 
interface h (Figure [2fc>) is chosen to be equal to the 
lattice elements crossing the interface between rein- 
forcement steel and concrete. The strains are related 
to the nominal stress a = (a n , cr s , a t ) T by the elasto- 
plastic stress-strain relationship 



a = D e (e - e c - s p ) 



(2) 
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Figure 2: Interface, (a) Tangential plane of interface with 
the local coordinate system n, s and t. (b) Cross-section of 
interface of thickeness h. 
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Figure 3: Yield surface: Mohr-Coulomb friction law com- 
bined with a cap. 

where D e is the elastic stiffness, e c = (e c , 0, 0) T is the 
eigenstrain describing the expansion of the corrosion 

product and e v = (e pn , s ps ,£ pt ) T is the plastic strain. 

The yield surface of the plasticity model consists 
of a Mohr-Coulomb friction law combined with an 
elliptical cap. The shape of the cap surface is adjusted 
so that a smooth transition between the two surfaces 
is obtained ( Figure [3]). This combin ation was initially 
proposed by |S wan and Seo (2000)| and further d evel- 
oped by |Dolarevic and Ibr ahimb egovic (2007)[ The 
yield function f depends on the normal stress a n and 

the shear stress norm a q = \J a 2 + a 2 as 

a q + aa n if a n o < cr n 



f={ 2 , K + fc 



P 2 



ifcr n <CTnO 



(3) 



where a is the frictional angle and f c is the compres 
sive strength. Furthermore, 

f3af c 



a = 



a(5 + y/1 + (5 2 a 2 



(4) 



where ft is the ratio of the short and long radii of the 
cap ellipse (Figure [3]). At the point where the two 
parts of the yield surface meet, the normal stress is 

rt 

. The rate of the plastic strains 



(J n () = , - 

in Equation [2] is 



A 



dg 
da 



(5) 



where g is the plastic potential and A is the plastic 
multiplier. In the present study, g is chosen to be very 
similar to the yield function /. The only difference is 
that a is replaced by the dilatancy angle i/j so that the 
magnitude of the normal plasticity strain generated 
during shear loading can be controlled. The plastic- 
ity model is completed by the loading and unloading 
conditions: 



/ < 0, A > 0, Xf = 



(6) 



This plasticity bond model is similar to the one devel- 
oped by |Lundgren (2005 a)[ However, in the present 
work, the response is assumed to be perfectly plastic, 
i.e. the shape of the yield surface is independent of 
the plastic strains. The bond model will be extended 
to hardening and softening in future work. The imple- 
mentation of the present model is simplified by intro- 
ducing a smooth transition between the cap and the 
frictional law. Thus, a special vertex stress return in 
the transition region is not required. 

2.2 Model for corrosion between concrete and rein- 
forcement 

The effect of corrosion is idealised as an eigenstrain 
e C9 which is determined from the free expansion of 
the corrosion product u c as e c = u c /h (Figure 3]). This 
expansion is related to the corrosion penetration depth 



^/r 2 + (2rx c -x 2 )(X c -l) - r (7) 

where A c is the ratio of the volume ratio of rust and 
steel. The corrosion penetration x c is related to the 
corrosion percentage p c as 



1-4/1- 



100 



(8) 



3 COMPARISON WITH EXPERIMENTAL RE- 
SULTS 

The lattice approach is used to model the experiments 
reported by |Lee et al. (2002)1 Th e geometry and load- 
ing setup of the experiments is shown in Figure E Re- 
inforcement bars (0 = 13 mm) embedded in concrete 
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Figure 4: Representation of the corrosion process as an ex- 
pansive layer of rust. 

cubes were initially subjected to corrosion and subse- 
quently pulled out. 
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Figure 5: Geometry and loading set-up for the by corrosion 
pull-out test reported by |Lee et al. (2002) [ The reinforce- 
ment bar with a diameter = 13 mm is placed eccentri- 
cally in ^/-direction in the concrete specimen. No lateral 
reinforcement is provided. 

The concrete used in the experiments is charac- 
terised by a Young's modulus of E c = 22.6 GPa, 
Poisson's ratio of u c = 0.17, a tensile strength of 
ft — 2.7 MPa and a compressive strength of f c = 
24.7 MPa. The Young's modulus of the reinforce- 
ment is E s = 183 GPa. In the present study ,the re- 
sponse of concrete, reinforcement and bond between 
concrete and reinforcement is modelled by the lattice 




Figure 6: Mesh for the lattice analysis. 

approach described earlier. The lattice for the analy- 
ses is shown in Figure [6l For the reinforcement and 
the interface between reinforcement and corrosion, 
the mesh is structured. For the concrete the lattice is 
random. 

Three analyses were performed. In the first anal- 
ysis the reinforcement was pulled out without initial 
corrosion. In the other two analyses, corrosion per- 
centages of p c = 3.2 and 16.8 % were considered 
before the pullout. Assuming uniform corrosion, the 
corrosion percentages were transformed according to 
Equation [8] to corrosion penetrations of x c = 105 and 
571 fim. With an expansion ratio of A c = 1.4, this 
gives, according to Equation^ free corrosion product 
expansions of u c = 41.6 and 214.8 /xm, respectively. 
With an interface element thickness of 1 mm, this 
gives eigenstrains of e c = 0.0416 and 0.2148. In all 
three analyses, the load F was controlled by the end 
slip in the form of relative horizontal displacements 
of nodes A and B as shown in Figure [5^. The results 
of the analyses are compared to the experimental re- 
sults in the form of average bond stress-slip curves 
shown in Figure [71 Here, the average bond stress was 
determined as r = Fj (tt0£) 9 where I = 60 is the em- 
bedded length. 

The pre-peak regime of the load- split curves ob- 
tained in the analyses is in very good agreement with 
the experiments. However, the post-peak responses 
obtained in the analyses deviate considerably from 
those reported in the literature. In particular, the anal- 
ysis of the corrosion free case exhibits a more brit- 
tle response than reported in the experiments. On the 
other hand, the load-slip curves with initial corrosion 
exhibit a more ductile response than reported in the 
experiments. This is surprising, since it is expected 
that cracking in the concrete cover should reduce the 
pressure at the interface and, thus, also the tangential 
stresses. Consequently, a more consistent pattern of 
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Figure 7: Comparison of predicted average bond stress-slip 
curves and experimental data reported by |Lee et al. (20 02) 
for three corrosion percentages p c = 0, 3.2 and 16.8 %. 

post-peak response might be expected across all three 
cases. More studies are required to explore this obser- 
vation. 

For the analyses without corrosion, the crack pat- 
terns for the peak bond stress and the maximum slip 
(presented in Figure [7]) are shown in Figure [H Crack 
patterns are visualised as those middle cross-sections 
of lattice elements in which the norm of the crack 
opening vector is greater than 10 /xm and increasing. 
Thus, only active cracks are presented. 

At the peak of the average bond stress-slip curve, 
the concrete cover is cracked at its thinnest section 
(Figured). With further slip, additional cracks initi- 
ate from the reinforcement and propagate radially into 
the specimen as shown in Figure [SJd. 

In Figure [9] the crack patterns are shown, for the 
two corrosion cases, at the end of the corrosion pro- 
cess. For both corrosion cases, cracking of the con- 
crete cover occurs before the pullout, which corre- 
sponds to the observations reported in the literature 
(ILee et al. 20021) . 

4 CONCLUSIONS 

In the present work a lattice approach is used to de- 
scribe the mechanical interaction of a corroding rein- 
forcement bar, the surrounding concrete and the inter- 
face between steel reinforcement and concrete. The 
cross-section of the ribbed reinforcement bar is taken 
to be circular, assuming that the interaction of the ribs 
of the deformed reinforcement bar and the surround- 
ing concrete is included in a cap-plasticity interface 
model. This lattice approach is capable of represent- 
ing many of the important characteristics of corrosion 
induced cracking and its influence on bond. The ide- 
alisation of the corrosion expansion as an eigenstrain 
allows for the modelling of corrosion induced crack- 
ing. Furthermore, the frictional bond law can model 




Figure 8: Crack patterns for the pullout analysis for the 
corrosion-free case at (a) peak and (b) end of the aver- 
age bond stress slip curve. Cracks initiate at the interface 
between reinforcement and concrete and propgate to the 
specimen surface. 

the decrease of the bond strength if the concrete is 
pre-cracked. Very good agreement with experimental 
results in the pre-peak regime of the bond stress-slip 
curves was obtained. More studies are required to in- 
vestigate the post-peak response of the bond stress- 
slip curves. Also, further studies will be performed to 
investigate the influence of the element length of the 
interface between reinforcement and concrete on the 
analyses results. Also, we will study the influence of 
the stiffness of the lattice elements on corrosion in- 
duced cracking and its interplay with lateral confine- 
ment. 
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